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ABSTRACT 

We extract synthetic photon spectra from first-principles particle-in-cell simulations of relativis- 
tic shocks propagating in unmagnetized pair plasmas. The two basic ingredients for the radiation, 
namely accelerated particles and magnetic fields, are produced self-consistently as part of the shock 
evolution. We use the method of Hededal & Nordlund (20051 and compute the photon spectrum 
via Fourier transform of the electric tar-field from a large number of particles, sampled directly from 
the simulation. We find that the spectrum from relativistic coUisionless shocks is entirely consistent 
with synchrotron radiation in the magnetic fields generated by Weibel instability. We can recover 
the so-called "jitter" regime only if we artificially reduce the strength of the electromagnetic fields, 
such that the wiggler parameter K = qBX/mc^ becomes much smaller than unity {B and A are the 
strength and scale of the magnetic turbulence, respectively). These findings may place constraints 
on the origin of non-thermal emission in astrophysics, especially for the interpretation of the hard 
(harder than synchrotron) low-frequency spectrum of Gamma-Ray Bursts. 

Subject headings: acceleration of particles — gamma rays: bursts — radiation mechanisms: non- 
thermal — shock waves 



1. INTRODUCTION 

Non-thermal photon spectra from Pulsar Wind Nebu- 
lae, jets from Active Galactic Nuclei, Gamma-Ray Bursts 
and Supernova Remnants are usually explained as syn- 
chrotron radiation from a power-law population of parti- 
cles, presumably accelerated in coUisionless shocks. The 
microphysical details of shock acceleration are still poorly 
known, however, and are the subject of active research. 

Particle-in-cell (PIC) simulations of colliding plasma 
shells have shown that W eibel instability ([Wsi bel 1959; 
Medvedev fc Loeb|[l999| iGruzinov fc Wa^an||1999| 
converts the free eiiergy of counter-stre aming fiows mto 
small scale ( skin-depth) magnetic fields (Nishikawa et al. 
^003) [20051 ISilva et al.|[2n03] [Frederiksen et al.||2004 



harder than expected from synchrotron radiation (e.g., 
Preece||1998|). Since PIC simulations self-consistently 



Hededal et al. 2004). The fields grow to sub-equipartition 
levels and defi ect and ra ndomize the bulk fiow, thus cre- 
ating a sho ck dSpit kovsky 2005, ,2008al [Chang et al.|2008 
Keshet et al. 2009). A few percent of the incoming par- 



ticles repeatedly scatter off the magnetic turbulence cre- 
ated by Weibel instability, and eventually populate a 
power-law high-energy tail in the particle spectrum be- 
hind the shoc k ( Spitkovsky 2008bl [Martins et al. [[2009b 
Sironi fc Spitko vsky 200g|T' 



provide both the strength and the spatial structure of 
electromagnetic fields, as well as the particle distribu- 
tion, it is possible to calculate the photon spectrum from 
first principles, thus determining whether the emission is 
more synchrotron-like or jitter-like. 

In this work, we present synthetic spectra extracted 
from PIC simulations of relativistic unmagnetized col- 
lisionless shocks. In §2 we show the simulated shock 
structure and particle energy spectrum. The numerical 
technique that we employ to compute the photon spec- 
trum is described in §3.1, and in §3.2 we show our results. 
In §4 we discuss the implications of our findings for the 
interpretation of the non-thermal radiation from astro- 
physical sources. 

2. SHOCK STRUCTURE AND PARTICLE ENERGY 
SPECTRUM 

We use the three-dimen sional (3D) elec t romagnetic 



PIC code TRISTAN-MP (Buneman 1993 Spitkovsky 



2005 1 to simulate a relativistic shock propagatmg mto an 



Since most of the magnetic power generated by Weibel 
instability is concentrated on scales as small as a few 
plasma skin depths, it has been speculated that the emis- 
sion mechanism in unmagnetized coUisionless shocks may 
be the so-called "jitter" radiation. Whereas the standard 
synchrotron emission applies to large-scale fields, the jit- 
ter regime is realized if the scale A of the turbulence is 
such that the w i ggler parameter K = qB X/mc^ <C 1 



( [Medvedev[[2000| |2006| |Fleishma!n l[2006a|b| . Jitter radi 



ation has been proposed as a solution for the so-called 
"line of death" puzzle in Gamma-Ray Burst (GRB) 
spectra, which below the peak frequency are sometimes 
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unmagnetized pair plasma. The shock is triggered by re- 
fiecting an incoming col d "upstream" fiow off a con duct- 
ing waU at a; = (e.g., iSironi fc Spitkovsk^[2009[ ) . The 
simulation is performed in the "wall" or "downstream" 
frame. The incoming fiow propagates along —x with 
Lorentz factor 70 = 15, and the shock moves along +x. 

To follow the shock evolution for longer times with 
fixed computational resources, we use a 2D simulation 
domain in the xy plane. In the case of an unmagne- 
tized 2D shock, only the in-plane components of the 
velocity, current and electric field, and only the out- 
of-plane component of the magnetic field are present. 
Each computational cell is initialized with 16 particles 
per species. The relativistic plasma skin depth for the 
upstream fiow (c/wp) is resolved with 10 cells and the 
simulation timestep is At = 0.045 Wp"^. The simulation 
box is lOOc/wp wide (along y) and, at the final Simula- 
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Fig. 1. — Shock structure and particle energy spectra at time 
u!pt = 2250. a) Number density in the simulation plane, normal- 
ized to the upstream density, b) Transversely-averaged density. 
c)-d) Magnetic and electric energy density in the simulation plane, 
normalized to the upstream kinetic energy density, e) Transversely- 
averaged magnetic (black) and electric (red) energy density, f) Lon- 
gitudinal phase space for positrons, shown as a 2D histogram (elec- 
tron phase space is nearly identical), g) Transverse phase space for 
positrons, h) Positron spectrum pdN/dp versus p (thick solid line; 
p = 7/3 = [(7/3x)^ + ilPy)^]^^^ is the particle 4-velocity) in the 
region x^h — 400c/a;p < x < + 400c/a;p (delimited by the ar- 
rows at the bottom of panel (g)), and the relative contributions 
of particles behind (thin solid line) and ahead (dotted line) of the 
shock. The dashed line corresponds to a power-law distribution 
dN/dp oc p~°' with a = 2.5. i) Positron spectrum pxdA^/dpx ver- 
sus px (with px = 7/3x) for particles with px > (red lines), and 
positron spectrum pydN/dpy versus py (with py = 7/3y) for par- 
ticles with Py > (blue lines). The line style coding is the same 
as in panel (h). The total spectrum around the shock (thick solid 
lines) is shifted upward by a factor of 10 for clarity. 



tion time ujpt = 4500, it is A500c/ujp long (along x). 

In Fig. [T] wc show the internal structure of the shock 
at time uipt = 2250. The 2D plots of number density 
(Fig. 1 1) , magnetic energy (Fig. [T];) and electric energy 
(Fig. Tl) show the filaments produced by Weibel insta- 
bility m the upstream region, with characteristic trans- 
verse scale ^ lOc/cUp. The magnetic filaments are ad- 
vected with the upstream fiow, so in the simulation frame 
their magnetic and electric components are comparable 
(black and red line in Fig. respectively). At the 
shock {xsh — 1030 c/wp at time ujpt = 2250), the fila- 
ments merge and the magnetic energy density peaks at 
~ 10% of the upstream kinetic energy density (Fig. Ilk) . 
The magnetic field decays farther downstream, where the 



field is confined within islands of typical scale ~ 20 c/wp 
(Fig#). 

The particle energy spectrum behind the shock (xsh — 
400 c/ujp < X < Xsh', thin solid line in Fig. [TJi) consists of 
a relativistic Maxwellian and a high-energy tail, which 
can be fitted as a power-law of index a — 2.5 (dashed 
line in Fig. Illi) with an exponential cutoff (Spitkovsky 
2008b I . In the upstream spectrum {xsh < x < Xsh + 
400 c/ujp; dotted line in Fig. [l|i), the unshocked beam 
populates the low-energy peak at p = 7/3 ~ 15, whereas 
the shock-accelerated returning particles with 7/3x > 
(see the hot, diS^use population in Fig. [if) appear as as 
a high-energy bump. Since at the highest energies the 
particles ahead of the shock account for nearly half of 
the total census (compare dotted and thin solid lines in 
Fig.fTri), the high-energy tail in the total spectrum (xgh — 
400 cfujp < X < a;sh+400 c/wp; thick solid line in Fig.[^) 
is significantly fiatter than in the downstream spectrum. 
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3. SYNTHETIC PHOTON SPECTRUM 
3.1. Numerical Technique 

the method introduced by IHededal 



summ arize 
and 
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2005[ ) (see also 
to extract syn- 



Marti ns et al. 2009a| 

thetic spectra from simulations of collisionless shocks 
The electric far-field from a particle wit h charge q, ve- 



locity V — (3c and acceleration = /3c is ( Jackson]|l999 | 



E{x,t) 



n X { (n — /3) X /3} 
(l-/3-n)3i? 



(1) 



where the unit vector n points toward the observer, at 
distance R from the emitting particle. Here, the quan- 
tity in square brackets is to be evaluated at the retarded 
time t' = t — R{t') / c. The photon spectrum is then com- 
puted via the Fourier transform of Poynting fiux associ- 
ated with the field in eq. ([TJ. The energy dW received 
per unit solid angle dfi (around the dire ction n) and p er 
unit frequency duj can be computed as ( Jackson]|1999 1 



dfidw A'k'^c 



n X {(n — (3) X [3} 
(l-/3-n)2 ' 



,iuj{t' -n-r{t')/c)^^l 



(2) 

where r{t') is the particle trajectory. Here, we neglected 
the Tsytovich-Razi n efi^ect due to the dispersi ve proper- 
tics of the plasma (Rybicki & Lightman' 1979). 

In PIC simulations we know the positions, velocities 
and accelerations of simulation particles with time reso- 
lution At = 0.045 ijjT^. In order to accurately compute 
the integral in cq. ([2| , we interpolate the orbit of the se- 
lected particles so that to achieve an effective timestep 
of 0.1 At. For a given choice of n, we can then integrate 
eq. ([2]) to obtain the photon spectrum from each parti- 
cle. Assuming that the far-fields by different particles are 
phase-uncorrelated, the total spectrum will be the sum 
of the spectra of individual particles. 

We have implemented eq. p| and tested it for the cases 
of synchrotron, bremsstrahlung and wiggler/undulator 
radiation, findin g good agreeme nt w ith analytic solu- 



tions . Follo wing |Hededal[ 
lund (2005), we have trace5~test particle 
netic turbulence with a prescribed wave spectrum. 



(|2005| and jHedcdal & Nord] 
latest particle orbits in mag- 



and 



we have verified the transition from the synchrotron to 
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Fig. 2. — Solid lines: photon spectrum from a power-law popu- 
lation of particles {dN/dp oc p~" between Pmia = 50 and pmax = 
2000, with a = 2.5) injected behind the shock (xgh — 400c/a;p < 
X < Xgjj). The particles are traced in the fields of the simulation 
at time ujpt = 2250. The photon spectrum from the same parti- 
cle population, but evolved in electromagnetic fields artificially re- 
duced by a factor of 10 (100, respectively), is shown as dot-dashed 
(dashed, respectively) lines. Red lines are for head-on emission 
(n = x), blue lines for edge-on emission (n = y). 

the jitter regime as the wiggler parameter K becomes 
smaller than unity. 

We now present the photon spectrum resulting from 
~ 10,000 particles moving in the electromagnetic fields 
self- consistently produced in our PIC simulations. We 
typically follow the particle trajectories over 3000 At = 
135 uj^^ , which is long enough to reach convergence in the 
shape of the spectrum, but short compared to the char- 
acteristic time of the shock evolution, such that the spec- 
trum we obtain may be regarded as "instantaneous." We 
compute the spectrum for both head-on emission (n = x, 
which we call "rix = 1" from now on) and edge-on emis- 
sion (n = y, "riy — 1" from now on). Our calculations 
are performed in the downstream fluid frame. An ad- 
ditional Lorentz transformation is required if the down- 
stream medium is moving with respect to the observer. 

3.2. Results 

We have studied the emission from relativistic shocks 
with two experiments. First, we have injected a power- 
law distribution of particles into the downstream region, 
and we have traced their orbits in a fixed snapshot of the 
fields of the simulation. The photon spectrum is then 
calculated with the technique described above. Second, 
we compute the spectrum from trajectories of particles 
extracted directly from the simulation. In this case, both 
the particle distribution and the electromagnetic turbu- 
lence are provided by the PIC simulation, and the result- 
ing emission will be the self- consistent spectrum from our 
relativistic shock. 

Fig. [2] shows the emission spectrum from a power- 
law population of particles behind the shock {xsh — 
400c/wp < X < Xsh)- We inject a 2D isotropic distri- 
bution with dN/dp oc p~°' {a = 2.5), between Pmin = 50 
and Pmax = 2000. The spectral index and the lower cutoff 
of the distribution are chosen to mimic the high-energy 
tail of the downstream particle spectrum in the simula- 



tion (thin solid line in Fig.[T|i). The injected particles are 
evolved in a fixed snapshot of the electromagnetic fields 
from the PIC simulation, at time Lo-pt — 2250. 

Since electric fields are negligible in the downstream 
medium (see Fig. [l]i-e), the resulting photon spectrum 
will probe the strength and structure of the magnetic 
fields, and it will clarify which regime - synchrotron or 
jitter - is appropriate to describe the particle emission. 
The solid lines in Fig. |2]show that the spectrum can be 
well approximated by two power-law segments. Regard- 
less of the observer's direction n (red line for = 1, 
blue line for riy = 1), the slope at the low frequencies 
is remarkably close to 2/3 (dotted line in Fig. |2|, as ex- 
pected for sync hrotron emissio n from a 2D particle dis- 



tribution (e.g., Jackson 1999). In N dimensions, the 
high-frequency slope should be —\a — (4 — N)]/2, which 
reduces to -(a-2)/2 = -0.25 for iV = 2 and a = 2.5, in 
agreement with our spectra (see dotted line). The sim- 
ilarity between the cases — 1 and Uy — 1 suggests, 
given the isotropy of the injected particle distribution, 
that the downstream magnetic fiuctuations are spatially 
isotropic, as seen in Fig. [ip.^ 

A transition to the jitter regime should appear when 
the wiggler parameter K becomes significantly smalle r 
than unity ( |Medved ev"2000l |2006| lFleishman|[2006a|b I . 
We have tried to artificially lower tne value of K by de- 
creasing the strength of the electromagnetic fields, by a 
factor of 10 (dot-dashed Hues in Fig. [2| and 100 (dashed 
lines in Fig.l2|. The high-frequency power-law decreases 
in intensity and shifts to lower frequencies, proportion- 
ally to the average magnetic field. The low- frequency 
spectrum becomes softer for decreasing B, approaching 
the fiat slope expected in the jitter regime when the 
shock is viewed edge-on. When we impose a magnetic 
field spectrum of the form B{k) (x S{k — feg) (here, fep 
is a fixed direction in fc-space) with wiggler parameter 
K <ti 1, we are able to r ecover the hard lo w-frequency 
slope (cx uj^) discussed by Medvedev (20001 for head-on 
emission from shocks. However, we do not observe it for 
the downstream turbulence self-consistently generated in 
the simulation, suggesting that the magnetic field fluctu- 
ations are not sufficiently ordered. 

In Fig.|3]we show the photon spectrum resulting from 
a sample of particles extracted directly from the PIC 
simulation, followed near cupt = 2250 in the time-varying 
electromagnetic flelds of the simulation. The selected 
particles start in the region Xsh — 400c/wp < a; < Xsh + 
400. The photon spectrum (thick solid lines in Fig.js] red 
for — 1, blue for Uy — 1) conflrms that the emission 
occurs in the synchrotron regime, as the 2/3 slope at the 
low frequencies suggests (black dashed line). 

Most of the low-frequency emission is powered by 
the thermal particles behind the shock, with 4-velocity 
p < 50 (dot-dashed lines in Fig. |3i). At high frequen- 
cies, the emission along y (blue s^id line in Fig. pa) is 
more powerful than along x (red solid line in FigTlSa). 
This reflects the fact that the highest energy par ticles 
are grazing the shock surface (Spitkovsky 2008b I and 



^ For a 3D distribution the low-frequency spectrum is oc oj^/"'. 

^ The slight difference at high frequencies between rix = 1 and 
riy = 1 is due to the residual electric fields at a; < (see Fig. Wi- 
e) , and it disappears if electric fields are neglected while computing 
the spectrum. 
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Fig. 3. — Photon spectrum (thick soUd Unes, in all panels) 
for particles extracted from the PIC simulation, evolved in time- 
varying electromagnetic fields around uipt = 2250. Red lines are for 
head-on emission (n = x), blue lines for edge-on emission {n = y). 
a) Dot-dashed lines only include the contribution of thermal par- 
ticles, with 4-velocity p < 50. b) Relative contribution of down- 
stream (thin solid lines) and upstream (dotted lines) particles, c) 
Time evolution of the photon spectrum, from uipt = 2250 (thick 
solid lines) to ujpt = 4500 (dot-dashed lines). Here, the case rix = 1 
has been shifted downward by a factor of 2 for clarity. 



they mostly contribute to the emission along y. The 
difference between Ux = f and riy = f in Fig. [3^ can 
be directly related to the difference between theparticle 
spectra px dA''/dpx and Py dN/dpy shown in Fig.fTl (thick 
solid lines, respectively red and blue). In fact, due to rel- 
ativistic beaming, an observer located along x will more 
likely detect the radiation from particles with Px/Py 3> 1 
(and px > 0), whereas the particles with Py/px ^ 1 (and 
Py > 0) will contribute more to the emission along y. 

As shown in Fig. [SJd, the spectrum of downstream 
particles (thin solid lines) dominates the total emission 
(thick solid lines) in the low-frequency bump, whereas 
the high frequencies are mostly powered by the spec- 
trum of upstream particles (dotted lines) . This is due to 
a combination of two effects: first, at the highest particle 
energies the upstream particles may dominate by num- 
ber (compare red dotted and thin solid lines in Fig. [TJ); 
second, the electromagnetic energy in the upstream re- 



gion is on average larger than in the downstream, since 
upstream electric fields are as strong as magnetic fields, 
and also because the field is built up ahead of the shock 
on a length scale larger than the decay scale in the down- 
stream (Fig.[T^). ft is worth pointing out that, at low fre- 
quencies, the spectra of both downstream and upstream 
particles can be fitted with the synchrotron slope 2/3 
(black dashed hues). 

Fig. ^ shows the time evolution of the total pho- 
ton spectrum, from ujpt = 2250 (thick solid lines) to 
tUpt = 4500 (dot-dashed hues). The upper cutoff in 
the particle energy spectrum grows linearly with time 



( |Spitkovsky|l2008b| |Sironi fc Spitkovsky|[2009| , and as a 
result the pnoton spectrum extends to higher frequen- 
cies. Also, the spectral intensity increases, especially at 
high frequencies, due to stronger electromagnetic fields 
and a larger number of accelerated particles. 

4. DISCUSSION 

We have computed synthetic photon spectra from 2D 
PIC simulations of relativistic coUisionless shocks by 
following a sample of simulation particles in the time- 
varying fields of the simulation. The low-frequency part 
of the spectrum scales as oc cj^/'^, as expected for syn- 
chrotron emission from a 2D particle distribution. Al- 
though the electromagnetic fluctuations generated by 
Weibel instability are on small (skin-depth) scales, the 
particle emission does not occur in the jitter regime. 

In retrospect, this is not surprising. The characteristic 
length scale of the magnetic turbulence is A > lOc/wp, 
and in the shock region, where most of the emission 
is produced, the magnetic energy reaches a fraction 
63 ~ 0.1 of the upstream bulk kinetic energy. This im- 
plies r^/{c/Lijp) = £3^^^ — 3, where is the relativistic 
Larmor radius of a particle moving with the upstream 
flow. It follows that K — X/{r^/-fo) — 3 70, so that 
the condition AT <C 1 for jitter radiation is unlikely to 
be satisfied, even for moderately relativistic shocks. In- 
deed, the photon spectrum we obtain from a 70 = 3 
shock is still consistent with synchrotron radiation. For 
electron-ion shocks, the wiggler parameter will be mp/m^ 
times larger, and we are even deeper in the synchrotron 
regime. Further downstream from the shock, although 
the magnetic field strength decays, the value of K does 
not significantly change, s ince short-wav elength modes 
are progressively damped ( Chang et al.| |2008j), and the 
characteristic scale of the turbulence increases. 

Although the results presented here apply to a 2D par- 
ticle distribution, the main conclusions should hold for 
3D configurations as well. There, we expect the low- 
frequency slope to be 1/3, as in the standard synchrotron 
radiation. If the GRB emission results from high-energy 
particles accelerated in relativistic unmagnetized shocks, 
it seems that resorting to the jitter radiati on is not a vi - 
able solution for the "line of death" puzzle ( Preece|1998 1 . 

At high frequencies, our results suggest that the contri- 
bution of upstream particles to the total emission, which 
is usually omitted in standard models, is not negligible. 
It causes the radiation spectrum to be flatter than the 
corresponding downstream spectrum, thus partly mask- 
i ng the contribution of downstream thermal particles 
(Giannios & Spitkovsky 2009). This could potentially 
explain the absence of clear signa tures of dow nstream 
thermal emission in GRB shocks (Band 19931. Simu- 



5 



lations extending to longer times (and higher particle 
energies) will help to clarify these issues. 

We remark that our calculations do not include ra- 
diative particle cooling, synchrotron self-absorption and 
inverse Compton radiation. Still, we have shown that 
the computation of synthetic spectra from self-consistent 
PIC simulations provides a powerful tool for studying 



the origin of astrophysical non-thermal emission. 

We thank J. Arons, A. Celotti, J. Kirk, S. Martins 
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Santa Barbara for hospitality. This research was sup- 
ported by NSF grants AST- 0807381 and PHY05-51164. 
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